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ABSTRACT 

Aims. We study the effects related to departures from non-local thermodynamic equilibrium (NLTE) and homogeneity in the atmo- 
spheres of red giant stars, to assess their influence on the formation of Ba II lines. We estimate the impact of these effects on the 
barium abundance determinations for 20 red giants in Galactic globular cluster NGC 6752. 

Methods. One-dimensional (ID) local thermodynamic equilibrium (LTE) and ID NLTE barium abundances were derived using clas- 
sical ID ATLAS9 stellar model atmospheres. The three-dimensional (3D) LTE abundances were obtained for 8 red giants on the lower 
RGB, by adjusting their ID LTE abundances using 3D~1D abundance corrections, i.e., the differences between the abundances ob- 
tained from the same spectral line using the 3D hydrodynamical and classical ID stellar model atmospheres. The 3D-1D abundance 
corrections were obtained in a strictly differential way using the 3D hydrodynamical and classical ID codes CO^'BOLD and LHD. Both 
codes utilized identical stellar atmospheric parameters, opacities, and equation of state. 

Results. The mean ID barium-to-iron abundance ratios derived for 20 giants are ([Ba/Fe]>iD lte = 0.24 ± 0.05(stat.) ± 0.08(sys.) and 
([Ba/Fe]>iD nlte = 0.05 ± 0.06(stat.) ± 0.08(sys.). The 3D-1D abundance correction obtained for 8 giants is small (~ +0.05 dex), 
thus leads to only minor adjustment when applied to the mean ID NLTE barium-to-iron abundance ratio for the 20 giants, 
([Ba/Fe]>3D+NLTE = 0. 10 ± 0.06(stat.) + 0.10(sys.). The intrinsic abundance spread between the individual cluster stars is small 
and can be explained in terms of uncertainties in the abundance determinations. 

Conclusions. Deviations from LTE play an important role in the formation of barium lines in the atmospheres of red giants studied 
here. The role of 3D hydrodynamical effects should not be dismissed either, even if the obtained 3D~1D abundance corrections are 
small. This result is a consequence of subtle fine-tuning of individual contributions from horizontal temperature fluctuations and dif- 
ferences between the average temperature profiles in the 3D and ID model atmospheres: owing to the comparable size and opposite 
sign, their contributions nearly cancel each other. This fine-tuning is characteristic of the particular set of atmospheric parameters and 
the element investigated, hence should not necessarily be a general property of spectral line formation in the atmospheres of red giant 
stars. 

Key words, stars: late-type - stars: abundances - stars: atmospheres - globular clusters: individual - techniques: spectroscopic 



.!^ 1. Introduction Although the GGC stars display a scatter in their Hght el- 

S^ ement abundances, there is generally no spread in the abun- 

;_( dances of iron-peak and heavier elements larger than the typ- 

.'^: Red giants in Galactic globular clusters (GGCs) carry a wealth ical measureme nt errors (^O.l dex). The only known ex cep- 

of important information about the chemical evolution of indi- tions are o) Cen (ISuntzeff& Kraflll 19961: iNorris et al.|[l99g) and 

vidual stars and their harboring populations. Owing to their in- M 54 (ICarrettaetalJ|2010b|), which do show noticeable star-to 



trinsic brightness, they are relatively easily accessible to high- star variations in the iron abundance. However, it is generally 
resolution spectroscopy, thus are particularly suitable for trac- accepted that they are not genuine GGCs but instead remnants 
ing the chemical evolution histories of intermediate age and of dwarf galaxies. The first cluster where significant start-to- 
old stellar populations. Unsurprisingly, a large amount of work star v ariation in heavy elemerit abundances was_ detected was 
has been done in t his direction in the past few decades (fo r M 15 (Snedenet al."199r, 2000; Otsuki et al. 2006; Sobeck et aiJ 
a review see, e.g., iGratton et al.l l2004t ICarretta at all l2010l) . |201 1). Roederer & Sneden (201 1) found that the abundances of 
which has resulted, for examp le, in the discoveries of abundance heavy elements La, Eu, and Ho in 1 9 red giants of M 92 indicate 
anti-correlations for Na-O dKraJ Ug^?; 'Gratton et al." "200?; that there are also significant star -to-star variations. The latter 
[Carretta et al. 2009a) Mg-Al (see, e.g., CaiTetta et al. 2009b), claim, however, was questioned by|Cohen] (201 1), who found no 
Li-Na (lPasauinietalJl2 005; Bonifaci oetani2007h. a nd correla- heavy element abundance spread larger than ~ 0.07 dex in 12 red 



tion for Li-O (iPasquini et al..,2005.: ,Shen et al.ll20'lOh . giants belonging to M 92. The primary formation channels of the 
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s-process elements are the low- and intermediate-mass asymp- 
totic giant branch (AGB) stars, thus the information about the 
variations in heavy element abundances may shed light on the 
importance of AGB stars to the chemical evolution of GGCs. 

Chemical inhomogeneities involving the light elements in 
GGCs are the result of the products of a previous generation 
of stars. The nature of the stars producing these elements, or 
'polluters' as they are often called, remains unclear. The main 
contenders are rapidly rota ting massive stars, th at pollute the 
cluster through their winds (iDecressin et al.ll2007l) or AGB stars 
(D'Ercole et al. 2011, and references therein). A second order 
issue is whether the "polluted" stars are coeval and only their 
photospheres are polluted or they are true second-generation 
stars formed from the polluted material. The evidence of multi- 
ple main-sequences and sub-giant branches in GGCs (see lPiotta 
12008 , 2009, for reviews) strongly supports the latter hypothesis, 
although some contamination of the photospheres may still be 
possible. 

Most of the abundance studies in GGCs have made the 
assumption of local thermodynamic equilibrium (LTE). Non- 
equilibrium effects may become especially important at low 
metallicity owing to th e lower opacities (e.g., overionization by 
UV photons; see, e.g., Asplund 2005; M ashonkina et alJl201 ll 
for more details). Deviations from LTE also occur because of 
the lower electron number density in the lower metallicity stel- 
lar atmospheres, which in turn decreases the electron collision 
rates with atoms and ions. Since most GGCs have metallicities 
that are significantly lower than solar, it is clearly desirable to 
derive abundances using the non-LTE (NLTE) approach. 

Nevertheless, real stars are neither stationary nor one- 
dimensional (ID), as assumed in the classical ID atmo- 
sphere models that are routinely used in stellar abun- 
dance work. A step beyond these limitations can be 
made by using three-dimensional (3D) hydrodynamical at- 
mosphere models that account for the three-dimensionality 
and non-stationarity of stellar atmospheres from first prin- 
ciples. Recent work has shown that significant differences 
may be expected between stellar abundances derived us- 
ing 3D hydrodynamical and classical ID model atmospheres 
(Collet etal. (2007, 2009); Gonzalez Hernandez et al. (200^ 
[Ramirez et al. (2009); Beharaetal. (2010); Dobrovolsk as et all 
(.2010.) ; .Ivanauskas et al, (.2010.) ; see also .Asplund (.20051) for a 
review of earlier work). These differences become larger at lower 
metallicities and at their extremes may reach 1 dex (!). 

It is thus timely to re-analyze in a systematical and homo- 
geneous way the abundances of various chemical elements in 
the GGCs, employing for this purpose state-of-the-art 3D hy- 
drodynamical atmosphere models together with NLTE analysis 
techniques. A step towards this was made in our previous work, 
where we derived ID NLTE abundances of Na, Mg, and Ba in 
the atmospheres of red giants belonging to GGCs MIO and M71 
("Mishenina e t al.ll2009l) . We found that in the case of the red gi- 
ant N30 in M71 the 3D-1D abundance corrections for Na, Mg, 
and Ba, were minor and did not exceed 0.02 dex. 

In this study, we extend our previous work and derive 
ID NLTE abundances of barium in the atmospheres of 20 red gi- 
ants that belong to the Galactic globular cluster NGC 6752. The 
analys is is done using the same techniques as in lMishenina et alJ 
(|2009|) . We also derive the 3D-1D LTE abundance corrections 
for the barium lines in 8 red giants and apply them to correct the 
ID barium abundances for the 3D effects. Finally, we quantify 
the influence of both NLTE and 3D-related efltects on the forma- 
tion of barium lines. 



The paper is organized as follows. In Sect. 2, we describe 
the observational material used in the abundance analysis. The 
procedure of barium abundance determinations is outlined in 
Sect. 3, where we also provide the details of the LTE/NLTE anal- 
ysis and the determination of the 3D-1D abundance corrections. 
A discussion of our derived results is presented in Sect. 4 and 
the conclusions are given in Sect. 5. 

2. Observational data 

We used reduced spectra of 20 red giants in NGC 6752 avail- 
able from the ESO Science ArchiveQ. The high resolution 
(R - 60000) spectral material was acquired with the UVES 
spectrograph at the VLT-UT2 (programme 65.L-0165(A), PI: 
F. Grundahl). Spectra obtained during the three individual ex- 
posures were co-added to achieve the average signal-to-noise 
ratio S /N ~ 130 at 600.0 nm. Observations were taken in the 
standard Die 346-1-580 nm setting that does not include the 
Ba II 455.403 nm resonance line. The other three Ba II lines 
at 585.369, 614.173, and 649.691 nm (see Table|2li are ah found 
in the upper CCD of the red arm covering the range 583-680 nm. 
More details of the spectra acquis ition and reduction procedure 
are provided bv lYong et al.l (l2005i) . All the red giants studied in 
this work are located at or below the red giant branch (RGB) 
bump. 



3. Abundance analysis 

3.1. Atmospheric parameters and iron abundances 

Continuum normalization of the observed spectra and equiva- 
lent width (EW) m easurements were m ade using the DECH28t1j 
software package ( Galazutd inov I fl992 ). where the EWs were 
determined using a Gaussian fit to the observed line profiles. 

Stellar model atmospheres used in the abundance de- 
terminations w ere ca lculated wi th the Linux port version 
(Sbordoneetal] l2004t !Sbordone1 l2005l) of the ATLAS9 code 
(Kurucz 1993 ), using the ODFNEW opacity distribution ta- 
bles from .Castelli & Kurucz (.2003.) . Models were computed 
using the mixing length parameter aMLT - 1-25 and micro- 
turbulence velocity of 1km s ', with the overshooting op- 
tion switched off. The LTE abundances were derived using the 
Linux port version (Sbordoneet al. 2004; Sbordone 2005) of 
the Kui-ucz W IDTHJ3 package (.Kurucz. 1993; .Kurucz . .2001 
ICastemll2005h . 

The effective temperature, Tgff, was determined under the 
assumption of excitation equilibrium, i.e., by requiring that the 
derived iron abundance should be independent of the excitation 
potential, x (Fig- [D upper panel). To obtain the value of surface 
gravity, log g, we required that the iron abundances determined 
from the Fe I and Fe II lines would be equal. The microturbu- 
lence velocity, ^t, was determined by requiring that Fe I lines of 
different EWs would provide the same iron abundance (Fig. [1] 
lower panel). The derived effective temperatures, gravities, and 
microturbulence velocities of individual stars agreed to within 
60 K, 0.2 dex, and . 16km s ', respectively, with those deter- 
mined bv lYongetan(l2005h . 

The LTE iron abundances for all stars in our sample were de- 
rived using 50-60 neutral iron lines (Table lATl note that the iron 
abundance derived from the ionized lines was requked to match 



http://archive.eso.org/eso/eso_archive_adp.html 

http://www.gazinur.com/DECH-software.html 

http://wwwuser.oat.ts.astro.it/castelli/sources/width9.html. 
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Table 1. Target stars, their adopted atmospheric parameters, and their derived iron abundances. 



Star 



RA 

(2000) 



Dec. 
(2000) 



Tcs, K 



dA/dx 
xlO-3, 
dex/eV 



^t, km s 



dA/dEW 
xlO-^ 
dex/pm 



iogg, 
[cgs] 



[Fe/H]i [Fe/H]„ 



A(Fe I) - A(Fe II), 
dex 



NGC 6752-1 


19:10:47 


-60:00:43 


4749 


-2.0 


1.37 


NGC 6752-2 


19:11:11 


-60:00:17 


4779 


-3.1 


1.35 


NGC 6752-3 


19:11:00 


-59:56:40 


4796 


6.9 


1.37 


NGC 6752-4 


19:11:33 


-60:00:02 


4806 


-4.9 


1.42 


NGC 6752-6 


19:10:34 


-59:59:55 


4804 


-7.6 


1.40 


NGC 6752-7 


19:10:57 


-60:00:41 


4829 


-11.1 


1.39 


NGC 6752-8 


19:10:45 


-59:58:18 


4910 


0.9 


1.31 


NGC 6752-9 


19:10:26 


-59:59:05 


4824 


5.8 


1.38 


NGC 6752-10 


19:11:18 


-59:59:42 


4836 


-1.4 


1.34 


NGC 6752-11 


19:10:50 


-60:02:25 


4870 


1.6 


1.33 


NGC 6752-12 


19:10:20 


-60:00:30 


4841 


3.5 


1.36 


NGC 6752-15 


19:10:49 


-60:01:55 


4850 


0.6 


1.35 


NGC 6752-16 


19:10:15 


-59:59:14 


4848 


-0.8 


1.35 


NGC 6752-19 


19:11:23 


-59:59:40 


4928 


^.0 


1.45 


NGC 6752-20 


19:10:36 


-59:56:08 


4929 


2.5 


1.25 


NGC 6752-21 


19:11:13 


-60:02:30 


4904 


-10.3 


1.34 


NGC 6752-23 


19:11:12 


-59:58:29 


4956 


1.7 


1.28 


NGC 6752-24 


19:10:44 


-59:59:41 


4948 


-2.8 


1. 19 


NGC 6752-29 


19:10:17 


-60:01:00 


4900 


1.9 


I.3I 


NGC 6752-30 


19:10:39 


-59:59:47 


4943 


-4.5 


1.18 



4.9 


1.95 


-1.53 


-1.50 


-0.03 


-0.2 


1.98 


-1.54 


-1.52 


-0.02 


-0.2 


2.03 


-1.59 


-1.60 


0.01 


0.8 


2.04 


-1.59 


-1.58 


-0.01 


0.3 


1.97 


-1.58 


-1.58 


0.00 


0.2 


2.10 


-1.77 


-1.75 


-0.02 


-0.3 


2.25 


-1.60 


-1.60 


0.00 


-0.9 


2.26 


-1.61 


-1.59 


-0.02 


0.3 


2.13 


-1.56 


-1.55 


-0.01 


0.9 


2.13 


-1.56 


-1.56 


0.00 


0.4 


2.15 


-1.60 


-1.63 


0.03 


-1.3 


2.19 


-1.57 


-1.56 


-0.01 


0.5 


2.06 


-1.63 


-1.63 


0.00 


0.1 


2.23 


-1.64 


-1.63 


-0.01 


-0.1 


2.33 


-1.56 


-1.59 


0.03 


0.0 


2.33 


-1.62 


-1.61 


-0.01 


1.9 


2.35 


-1.57 


-1.55 


-0.02 


1.3 


2.28 


-1.63 


-1.63 


0.00 


-2.0 


2.24 


-1.68 


-1.63 


-0.05 


-0.3 


2.42 


-1.64 


-1.62 


-0.02 



mean 



-1.60 
0.05 



-1.60 
0.05 



Note: Fe II abundances were adjusted to match the abundances determined from Fe I lines, in order to estimate the surface gravities of the target 
stars. The difference between the corresponding abundance ratios, [Fe/H];, and [Fe/HJu, is thus only indicative of the goodness of the gravity 
estimates. 



that of neutral iron, i.e., to obtain the estimate of surface grav- 
ity, thus it is not an independent iron abundance measurement). 
To minimize the impact of NLTE effects on the iron abundance 
determinations, we avoided neutral iron lines with the excita- 
tion potential ;if < 2.0 eV. Oscillator strengths and damping con- 
stants for all iron lines were retrieved from the VALD database 
(|Ktipka e t al. 2000). The obtained iron abundances are provided 
in Table[T] The contents of the table are as follows: the star iden- 
tification and its coordinates are given in Cols. 1-3, effective 
temperatures and iron abundance deiivatives relative to the exci- 
tation potential are in columns 4 and 5, respectively, the adopted 
microturbulence velocity and iron abundance derivative relative 
to the equivalent width are in columns 6 and 7, respectively, the 
adopted values of log^ are in column 8, iron abundances ob- 
tained from Fe I and Fe II lines are in columns 9 and 10, respec- 
tively, and the difference between them is in column 1 1 . 

The mean iron abundance obtained for the 20 stars is 
([Fe/H]) = -1.60 + 0.05, which is in excellen t agree ment with 
[Fe/H] = -1.62 + 0.02 obtained bv lVong et al.l (l2()05h . 

3.2. One-dimensional LIE abundances of barium 

One-dimensional (ID) LTE barium abundances were derived 
from the three Ba II lines centered at 585.3688 nm, 614. 1730 nm, 
and 649.6910nm. Damping constants and other atomic param- 
eters of the three barium lines are provided in Table The line 
equivalent widths were measured with the DECH20T software 
(Table[3] columns 2-A). Hyperfine splitting of the 649.6910nm 
line was not taken into account in the ID LTE analysis. The de- 
rived barium abundances and barium-to-iron abundance ratios 
are given in Table [3] columns 5 and 7, respectively. 

We note that the barium line at 614. 1730nm is blended with 
a neutral iron line located at 6 14. 171 3 nm. To estimate how this 



Table 2. Atomic parameters of the barium lines used in this 
work. 



, nm 



X\ eV log gf^ log Jra 



log K 



los -^'' 



:vfl_ 



585.3688 
614.1730 
649.6910 



0.604 
0.704 
0.604 



-1.000 
-0.076 
-0.377 



8.20 
8.20 
8.10 



-5.460 
-5.480 
-5.480 



-7.190 
-7.470 
-7.470 



Wiese & Martin ("1980^; '' na tural broadening constant, from 
Mashonkina & Bikma^ J1996I) : '^ Stark broadening constant, from 
Kuokaetal. 12000); "^ van der Waals broadening constant, from 
iKorotinet ah l i2011i) 



affects the accuracy of the abundance determination, we synthe- 
sized the barium 614.1730nm line with and without the blend- 
ing iron line, for all stars in our sample. The comparison of 
the equivalent widths of blended and non-blended lines reveals 
that the contribution of the iron blend never exceeds ~ 2.4%, 
or < 0.05 dex in terms of the barium abundance. The contri- 
bution of the iron blend to the EW of the 614.1730nm line 
was thus taken into account by reducing the measured equiva- 
lent widths of this barium line by 2.4 % for all stars. We would 
like to point out, however, that in the ID NLTE analysis the bar- 
ium abundances were derived by fitting the synthetic spectrum 
to the observed line profile, thus the influence of the iron blend 
at 614.1713nm was properly taken into account. 

Assessment of the abundance sensitivity on the atmospheric 
parameters yields the following results: 

- Change in the effective temperature by ±80 K leads to a 
change in the barium abundance measured from the three 
Ba II lines by +0.03 dex; 
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Table 3. Measured equivalent widths of the barium lines and the derived barium abundances in individual red giants in NGC 6752. 



Star 



EW (585.4nm), 
pm 



£iy(614.2nm), 
pm 



EW (649.7nm), 
pm 



A(Ba)iD LTE, 

dex 



A(Ba)iD NLTE, 

dex 



[Ba/Fe]iDLTE, 
dex 



[Ba/Fe]iDNLTE, 
dex 



6752-1 


8.21 


12.54 


12.87 


0.92 


0.63 


0.28 


-0.01 


6752-2 


8.43 


12.45 


12.08 


0.91 


0.65 


0.28 


0.02 


6752-3 


7.21 


11.84 


11.65 


0.75 


0.50 


0.17 


-0.08 


6752-4 


8.25 


12.07 


12.26 


0.86 


0.68 


0.28 


0.10 


6752-6 


8.23 


11.96 


12.80 


0.87 


0.60 


0.28 


0.01 


6752-7 


7.13 


10.72 


11.32 


0.60 


0.44 


0.20 


0.04 


6752-8 


7.01 


11.29 


11.37 


0.84 


0.72 


0.27 


0.15 


6752-9 


7.50 


11.65 


12.02 


0.85 


0.70 


0.29 


0.14 


6752-10 


7.57 


11.58 


12.16 


0.87 


0.67 


0.26 


0.06 


6752-11 


7.59 


11.62 


11.13 


0.82 


0.61 


0.21 


0.00 


6752-12 


7.01 


11.44 


11.34 


0.76 


0.61 


0.19 


0.04 


6752-15 


7.13 


11.24 


11.35 


0.78 


0.60 


0.18 


0.00 


6752-16 


7.36 


11.51 


11.59 


0.80 


0.58 


0.26 


0.04 


6752-19 


6.84 


10.76 


11.03 


0.69 


0.57 


0.16 


0.04 


6752-20 


6.77 


11.16 


10.81 


0.87 


0.75 


0.26 


0.14 


6752-21 


7.29 


10.74 


11.09 


0.82 


0.64 


0.27 


0.09 


6752-23 


7.17 


10.71 


11.21 


0.88 


0.72 


0.28 


0.12 


6752-24 


6.03 


9.79 


9.57 


0.68 


0.53 


0.14 


-0.01 


6752-29 


6.35 


10.06 


10.52 


0.67 


0.51 


0.18 


0.02 


6752-30 


6.39 


9.92 


10.57 


0.82 


0.66 


0.29 


0.13 



mean 
a 



0.80 
0.09 



0.62 
0.08 



0.24 
0.05 



0.05 
0.06 



-1.2 

-1.4 

ffi -1.6 

El 
-1.8 

-2.0 



-1,4 
ffi -1.6 
^ -1.8 

-2.0 



1 ' 1 
■ NGC 6752-23 


' 1 


' 1 


o 
1 1 1 


d[Fe/H]/dx=1.7 
1 1 


xlo'dex/eV . 
1 1 



X,eV 



--% 


1 1 1 1 
o 


1 

OT) 


V-O^^ 


' 1 

©o 


•^. 


1 1 1 1 
I 


. . . \ . . . 


1 


oOq * • 

d[Fe/H]/d£PF= 
, . , 1 , , , 1 


• 
1.9x10 


• 

dex/pm 
1 , , , 



6 
EW, pm 



10 



12 



Fig. 1. [Fe/H] abundance ratios derived from Fe I lines for the 
star NGC 6752-23, plotted versus the excitation potential (top) 
and line equivalent width (bottom). Linear fits to the data are 
shown as solid lines. Note that there is a slight correlation be- 
tween the best-fit slopes in the two panels: adjusting the temper- 
ature by zeroing the slope in the upper panel slightly changes 
the slope in the lower panel, while changes to the microturbu- 
lence velocity influence the slope in the upper panel. This cor- 
relation does not allow us to obtain zero-valued slopes in both 
panels simultaneously. The adopted atmospheric parameters of 
this star are T^ff = 4956 K, logg = 2.35, ^t - 1.28km s"\ and 
[Fe/H] = -1.57. 



- Change in the surface gravity by ±0. 1 dex changes the bar- 
ium abundance by +0.02 dex; 



- Change in the micro turbulence velocity, ^t, by 
changes the barium abundance by +0.07 dex. 

Since barium lines in the target stars are strong and situated 
in the saturated part of the curve of growth, it is unsurprisi ng that 
the uncertainty in the microturbulence velocity is the largest con- 
tributor to the uncertainty in the derived barium abundance. The 
total contribution from the individual uncertainties in Teff , log g, 
and ^t leads to the systematic uncertainty in the barium abun- 
dance determinations of ~ 0.08 dex. We note, however, that the 
latter number does not account for the uncertainty in the equiva- 
lent width determination and thus only provides the lower limit 
to the systematic uncertainty (e.g., 5 percent in the equivalent 
width determination leads to the barium abundance uncertainty 
of -0.1 dex). 

The obtained mean ID LTE barium abundance for the sam- 
ple of 20 stars in NGC 6752 is <A(Ba)>iD lte = 0.80+0.09+0.08 
and the barium-to-iron ratio is ([Ba/Fe])io lte - 0-24 + 0.05 ± 
0.08. In both cases, the first error is a square root of the variance 
calculated for the ensemble of individual abundance estimates of 
20 stars. The second error is the systematic uncertainty in the at- 
mospheric parameter determination. The difference between the 
individual barium abundances derived in a given star using the 
three barium lines is always below ~ 0.1 dex. 

3.3. One-dimensional NLTE abundances of barium 

The one dimensional (ID) NLTE abundances of barium were 
determined using the version of the ID NLTE spectral synthesis 
code MULTI (Carlsson 1986) modified by Korotin et al. ( 1999|). 
The model atom of barium u sed in the NLTE spectra l synthesis 
calculations was taken from lAndrievskv et al.l (1200911 To sum- 
marize briefly, it consisted of 3 1 levels of Ba I, 101 levels of Ba II 
(n < 50) and the ground level of Ba III. In total, 9 1 bound-bound 
transitions were taken into account between the first 28 levels of 
Ball (n < 12, Z < 5). Fine structure was taken into account for 
the levels 5d^D and 6p^P", according to the prescription given 
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Fig. 2. Fit of the synthetic NLTE profiles of the Ba II lines (solid line) to the observed spectrum of the red giant NGC6752-23 (dots). 
Synthetic line profiles corresponding to the abundances 0.1 dex higher/lower than the best-fit value are shown as dashed lines. 



in lAndrievskv et all (l2009h . We also accounted for the hyperfine 
splitting of the barium 649. 6910 nm line. Isotopic splitting of the 
barium lines was not taken into account. Owing to the low ion- 
ization potential of neutral barium (~ 5.2 eV), Ba II is the dom- 
inant ionization stage in the line-forming regions of investigated 
stars, with n(Ba I)/n(Ba II) < lO"'*. It is therefore safe to assume 
that none of the Ba I tra nsitions may noticeably change the level 
populations of Ba II (cf. lMashonkina et aljl999 ). Further details 
about the barium model atom, the assumptions used , and im pli- 
cations involved can be found in lAndrievskv et al.l (l2009l) and 
iKorotinet £11(1201 lb . 

The solar abundances of iron and barium were assumed to 
be log A(Fe)0 = 7.50 and log A(Ba)0 = 2.17 respectively, on the 
scale where log A(H)0 = 12. These abu ndances w ere determined 
using the Kurucz Solar Flux Atlas ( Kurucz et alj|l984l) and the 
same NLTE approach as applied in this study. 

A typical fit of the synthetic line profiles to the observed 
spectrum is shown in Fig. |2] where we plot synthetic and ob- 
served profiles of all three barium lines used in the analysis. The 
elemental abundances and barium-to-iron abundance ratios de- 
rived for the individual cluster giants are provided in Table [3] 
(columns 6 and 8, respectively). 

The mean derived ID NLTE barium- to-iron ratio for the 20 
cluster red giants is <[Ba/Fe])iD nlte = 0.05 + 0.06 + 0.08. The 
first error is the square root of the variance in ([Ba/Fe])iQ ^lte 
estimates obtained for the ensemble of 20 stars, thus measures 
the star-to-star variation in the barium-to-iron ratio. The second 
error is the systematic uncertainty resulting from the stellar pa- 
rameter determination (see Section 3.1). The individual line-to- 
line barium abundance scatter was always significantly smaller 
than 0. 1 dex. 

We find that barium lines generally appear stronger in NLTE 
than in LTE, which leads to lower NLTE barium abundances. 
This is in accord with the results obtained bv lShort & Hauschildtl 
(120061) for the metallicity of NGC 6752, and similar to th e 
trends obtained for cool dwarfs by iMashonkina et al.l (Il999l) . 
The NLTE-LTE abundance corrections for the three individual 
barium lines are always very similar, with the differences being 
within a few hundredths of a dex. 



Table 4. Parameters of the 3D hydrodynamical CO^BOLD atmo- 
sphere models used in this work. 
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3.4. 3D-1D barium abundance corrections 



We h ave used the CO^BOLD 3D hydrodynamical (Frevta g et al.l 
12012!) and LHD ID hydrostatic ( Caff^au & Ludwig ..2007.) stellar 
atmosphere models to investigate how strongly the formation of 
barium lines may be affected by convective motions in the stellar 
atmosphere. The CO^BOLD code solves the 3D equations of radi- 
ation hydrodynamics under the assumption of LTE. The model 
assumes a cartesian coordinate grid. For a detailed description of 
the CO^BOLD code and its applications, we refer to lFrevtaget al.l 
(120120 . 

Since we did not have CO^BOLD models available for the en- 
tire atmospheric parameter range covered by the red giants in 
NGC 6752, we estimated the importance of 3D hydrodynami- 
cal effects only for stars on the lower RGB. For this purpose, 
we used a set of 3D hydrodynamical CO^BOLD models with 
Teif = 5000 K and logg = 2.5, at four different metallicities, 
[M/H]= 0.0, -1.0, -2.0, and -3.cB Allowing for the eiTor mar- 
gins of ~ lOOK in the effective temperature and ~ 0.25 dex in 
gravity, we assumed that the effective temperature and gravity 
of this model set is representative of the atmospheric parameters 
of the stars NGC 6752-08, and NGC 6752-19 to NGC6752-30 
(8 objects, see Table [TJ. For these stars, the extreme deviations 
from the parameters of the 3D model are AT^eff ~ llOK and 
A logg ~ 0.26. These differences would only have a marginal 
effect on the uncertainty in the abundance estimates, i.e., the 
systematic uncertainty for the 3D barium abundance deriva- 
tions would only increase from +0.08 dex quoted in Sect. 13.21 
to +0.10 dex. 

The 3D hydrodynamica lmodels were taken from the CIFIST 
3D model atmosphere grid (iLudwig et al.ll2009h . The model pa- 



"^ The models at four metallicities were needed to interpolate the 3D- 
ID abundance corrections at the metallicity of NGC 6752. 
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Fig. 3. Top panel: Temperature stratification in a single snap- 
shot of the 3D hydrodynamical CO^BOLD model at the metallicity 
[M/H]= -2. Gray shaded area shows the temperature probabil- 
ity density on a logarithmic scale, with darker shades meaning a 
higher probability of finding a particular temperature value in the 
3D model simulation box. The solid red line shows mean tem- 
perature stratification of the 3D model and the dashed red line 
is the ID LHD model temperature stratification. Horizontal bars 
show the optical depth intervals where 90% of the line equiv- 
alent width is formed: black bars correspond to the 3D model 
while blue dashed correspond to the ID. Numbers next to the 
bars designate the wavelength of the given Ba II line in nm. 
Bottom panel: RMS value of horizontal temperature fluctuations 
in the 3D model (black line) and temperature diff'erence between 
the mean 3D and ID models (blue dashed line). 



rameters are summarized in Tabled The physical size of the 3D 
model box was chosen so that it would accommodate at least ten 
convective cells in the horizontal plane. Monochromatic opaci- 
ties used in the model calculations were grouped into five opac- 
ity bins for [M/H] = 0.0 and six opacity bins for [M/H] - 
-1.0,-2.0,-3.0 models. The ID LHD models were calculated 
for the same set of atmospheric parameters using the same equa- 
tion of state, opacities, and chemical composition as in the case 
of the 3D hydrodynamical models. 

To illustrate the differences between the 3D hydrodynami- 
cal and ID classical stellar model atmospheres, we show their 
temperature stratifications at the metallicity of [M/H] - -2.0, 
which is the closest to that of NGC 6752 (Fig. [3] upper panel). 
In the same figure, we also indicate the typical formation depths 
of the three barium lines. It is obvious that at these depths, the 
temperature of the 3D hydrodynamical model fluctuates very 
strongly, especially in the outer atmosphere, as indicated by the 
RMS value of horizontal temperature fluctuations (AFrms - 
-^((r - ro)^).i-,y.r, where To is the temporal and horizontal tem- 
perature average obtained on surfaces of equal optical depth). 
As we see below, differences in the atmospheric structures lead 
to differences in the line formation properties and henceforth to 
differences in barium abundances obtained with the 3D hydro- 
dynamical and ID classical model atmospheres. 



Twenty 3D snapshots (i.e., 3D model structures at different 
instants in time) were selected to calculate the Ba II line pro- 
files. The snapshots were chosen in such a way that the statistical 
properties of the snapshot sample (average effective temperature 
and its rm.s value, mean velocity at the optical depth unity, etc.) 
would match as close as possible those of the entire ensemble of 
the 3D model run. The 3D line spectral synthesis was performed 
for each individual snapshot and the resulting line profiles were 
averaged to yield the final 3D spectral line profile. 

The influence of convection on the spectral line formation 
was estimated by means of 3D-1D abundance corrections. The 
3D-1D abundance correction is defined as the difference be- 
tween the abundance A{Y) derived for a given element Y from 
the same observed spectral line using the 3D hydrodynamical 
and classical ID model atmospheres, i.e., A3D-1D = MY)3d - 
A(Y)iu (Caffauetal. 2011). This abundance correction can be 
separated into two constituents: (a) correction owing to the hori- 
zontal temperature inhomogeneities in the 3D model, A3d-(3d> = 
MY)3D-A{Y)(3D), and (b) correction owing to the differences be- 
tween the temperature profiles of the average (3D) and ID mod- 
els, A<3D)-iD = A{Y)(3D) -MY)iD- Abundances corresponding to 
the subscript (3D) were derived using the average 3D models, 
which were obtained by horizontally averaging 3D model snap- 
shots on surfaces of equal optical depth. Spectral line profiles 
were calculated for each average (3D) structure corresponding 
to individual 3D model snapshots. These line profiles were aver- 
aged to yield the final (3D) profile, which was used to derive the 
A<3D>-iD abundance corrections. The full abundance correction 
was then A3D-1D = A3d-(3d> + A<3d)-id- Spectral line synthesis 
for all three models, i.e., 3D, (3D), and ID, was made using the 
Lin£or3Dcodfl 

The barium lines in the target stars are strong (cf. Table |3) 
and thus the derived 3D-1D abundance corrections are sensitive 
to the microturbulence velocity, ^i, of the comparison ID model. 
The 3D-1D abundance corrections were therefore calculated us- 
ing the equivalent widths and microturbulence velocities of the 
target stars derived in Sect. 13. H and 1X2] Furthermore, cubic inter- 
polation was used to interpolate between the 3D-1D abundance 
corrections derived at four different metallicities to obtain its 
value at the metallicity of the cluster, [Fe/H] - -1.6. The cubic 
interpolation between the four values of metallicities was chosen 
because of the nonlinear dependence of the 3D-1D abundance 
corrections on metallicity. The results are provided in Table |5] 
which contains the A3D-1D and A3d-<3d> abundance corrections 
for the three individual Ba II lines (columns 2-4), the 3D-1D 
abundance correction for each star (i.e., averaged over the three 
barium lines, column 5), the microturbulence velocity used with 
the ID comparison model (column 6, from Sect.[T]l, the 3D LTE 
barium abundances (column 7), the 3D LTE barium-to-iron ratio 
(column 8), and finally both the ID NLTE barium-to-iron ratio 
before (column 9, from Sect. I3.3l l and after correction for the 3D 
effects (column 10). 

Abundance corrections are sensitive to the choice of the ID 
microturbulence velocity and line strength, therefore stars with 
very similar atmospheric parameters may have different abun- 
dance corrections. This is, for example, the case for NGC 6752- 
19 and NGC 6752-30. These two stars have largest and smallest 
microturbulence velocities in the entire sample, respectively, and 
NGC6752-19 has sHghtly stronger barium lines than NGC 6752- 
30 (Table |3]l. This leads to noticeably different abundance cor- 
rections, despite both stars having very similar effective temper- 
atures and gravities (Table|5). 



http://www.aip. de/"'mst/Linfor3D/linfor_3D_manual.pdf 
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Table 5. The 3D LTE barium abundances of red giants in NGC 6752. 



Star 



Aline 



""'= , dex 



BaII5854 



A' 

D' "3D-<3D> 

BaII6142 



BaII6497 



A3D-1D 
dex 



km s"' 



A(Ba) 
3D LTE 



[Ba/Fe] 
3D LTE 



[Ba/Fe] 
ID NLTE 



[Ba/Fe] 
3D+NLTE 



6752-08 


0.05 0.00 


0.08 


-0.01 


0.06 


-0.03 


0.06 


1.31 


0.90 


0.33 


0.15 


0.21 


6752-19 


0.13 0.08 


0.18 


0.11 


0.15 


0.07 


0.15 


1.45 


0.84 


0.31 


0.04 


0.19 


6752-20 


0.01 -0.04 


0.04 


-0.05 


0.02 


-0.07 


0.02 


1.25 


0.89 


0.28 


0.14 


0.16 


6752-21 


0.07 0.02 


0.10 


0.02 


0.08 


-0.01 


0.08 


1.34 


0.90 


0.35 


0.09 


0.17 


6752-23 


0.02 -0.03 


0.06 


-0.02 


0.04 


-0.05 


0.04 


1.28 


0.92 


0.32 


0.12 


0.16 


6752-24 


-0.02 -0.07 


-0.01 


-0.08 


-0.03 


-0.10 


-0.02 


1.19 


0.66 


0.12 


-0.01 


-0.03 


6752-29 


0.05 0.00 


0.08 


0.01 


0.06 


-0.02 


0.06 


1.31 


0.73 


0.24 


0.02 


0.08 


6752-30 


-0.03 -0.08 


-0.01 


-0.09 


-0.02 


-0.11 


-0.02 


1.18 


0.80 


0.27 


0.13 


0.11 



mean: 
sigma: 



0.83 
0.09 



0.28 
0.07 



0.08 
0.06 



0.13 
0.08 



4. Discussion 

4.1. One-dimensional LTE and NLTE barium abundances in 
NGC 6752 

Generally, the mean ID LTE barium-to-iron abundance ratio ob- 
tained in this work, ([Ba/Fe])iD lte - 0.24 + 0.05 (random) 
+ 0.08 (systematic), agrees well with the ID LTE abundance 
ratios derived for this cluster by other authors. For example, 
[James et al. (2004) derived ([Ba/Fe]) = 0.1 8+0. 11 from 9 main- 
sequence and 9 subgiant star s in NGC 6752. We note that the 
mean barium-to-iron ratio of [James et all (l2004b based on the 
measurements of only subgiant stars is ([Ba/Fe]) = 0.25 + 0.08, 
i.e., in this case the agreement with our LTE estimate is even 
better. The barium-to-i ron ratios obtained by Norris & D a Costal 
(IT99I and lYonget all (12005) are lower, ([Ba/Fe]) = 0.00 + 0.13 
and ([Ba/Fe]) = -0.06 + 0.13, respectively. 

The disagreement between the results of [Yong et al.l ([20051) 
and those obtained here is somewhat concerning, especially 
since both studies were based on the same set of UVES spec- 
tra, while the atmospheric parameter s and iron abundan ces of 
individual stars employed by us and [Yong et al] ( [2005h agree 
very well (Sect. 13. 11 1. Moreover, the compari son of the equiva- 
lent width measurements obtained by us and [Yong et aP ([2005h 
also shows good agreement. One would thus also expect good 
agreement in the derived barium abundances - which is unfor- 
tunately not the case. We therefore felt it was important to look 
into the possible causes of this discrepancy. 

To this end, we first obtained the ID LTE barium abundances 
using the MULTI code. This independent abundance estimate was 
made using the same procedure as for the ID NLTE abundance 
derivations, i.e., by fitting the observed and synthetic line pro- 
files of the three Ba II lines, with the difference that in this case 
the line profile calculations performed with MULTI were done 
under the assumption of LTE. The mean barium-to-iron abun- 
dance ratio obtained in this way, ([Ba/Fe]) - 0.22 + 0.06 + 0.08, 
agrees well with the value derived in Section 3.2 (([Ba/Fe]) - 
0.24 ± 0.05 + 0.08). 

In their abundance determinations, Yong et aP ([20051) used 
an older version of the ATLAS models ([Kuruczll993l) . The differ- 
ences between these ATLAS models and those used in our work 
is that (a) differe nt opacity tables were use d in the two cases (i.e., 
ODFNEW from [Castell i & Kur uczl[2003L with our models), and 
(b) the ATLAS models of Kumcz jl993h were calculated with the 
overshooting parameter switched on, while in our case the over- 
shooting was switched off. To check the influence of these dif- 
ferences on the abundance derivations, we obtained the ID LTE 
barium abundance using the older atmosphere models of [Kuruca 



([1993[) . with the atmospheric parameters and iron abundances 
derived in Sect. 13.1! In this case, the mean derived barium-to- 
iron abundance ratio was ([Ba/Fe])iD lte = 0.23 + 0.05 ± 0.08, 
i.e., the effect of differences in the model atmospheres was only 
~ 0.01 dex. The change in the barium abundances owing to dif- 
ferences in the atomic parameters (line broadening constants, 
oscillator strengths) used in the two studies was more signifi- 
cant, i.e ., the abundances d erived by us using the atomic param- 
eters of [Yong et al.l ([2005[) were -0.1 dex lower. However, this 
still leaves a rather large discrepancy, ~ 0.15 dex, between the 
barium-to-iron ratios obtained by us and Yong et al. (2005), for 
which we unfortunately cannot find a plausible explanation. 

As in the case of the ID LTE abundances, the extent of the 
star-to-star variations in the derived ID NLTE barium-to-iron ra- 
tio, ([Ba/Fe])[DNLTE = 0.05 ± 0.06 + 0.08, is small and can be 
fully explained by the uncertainties in the abundance determina- 
tion. The ID NLTE barium-to-iron ratio derived here is similar 
to the value ([Ba/Fe]), nNrxF. - 0.09 + 0.20 obtained for two 
red giants in MIO by Mishenina et al. (2009). The elemental ra- 
tios obtained in the two studies are thus very similar, although 
one sh ould keep in mind that the estimate of Mishenina et al] 
([2009[) is based on only two stars. The metallicities of the two 
clusters are very similar too, [Fe/H] - -1 .56 in the case of MIO 
(Harris 1996, 2010 edition) and [Fe/H] = -1.60 for NGC 6752 
(Sect. 13.1b . Galactic field stars typically show no pronounced 
dependence of [Ba/Fe] on metallicity, although the scatter at any 
given metallicity is large ([Sneden et al J [20081) . One may there- 
fore conclude that, taken into account the high [Ba/Fe] spread 
in field stars, the ([Ba/Fe]) ratio derived here is comparable to 
those seen in Galactic field stars and other globular clusters of 
similar metallicity. 

4.2. The 3D-corrected 1D NLTE barium abundance in 
NGC 6752 

The 3D-1D barium abundance corrections obtained for the eight 
stars in NGC 6752 (see Section 3.4 above) provide a hint of 
the net extent to which the 3D hydrodynamical effects may 
influence spectral line formation (and thus, the abundance de- 
terminations) in their atmospheres (Table |5)- In the case of 
all red giants investigated, the corrections are small, -0.03 to 
+0. 15 dex, and the mean abundance correction for the eight stars 
is A3D-1D = 0.05. We note though that the individual contri- 
butions to the abundance correction, A3d-(3d> and A(3d>-id, are 
substantial (~ ±0.1 dex) but often because of their opposite sign 
nearly cancel and thus the resulting abundance correction is sig- 
nificantly smaller (Table|5]l. This clearly indicates that the role of 
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convection-related effects on the spectral line formation in these 
red giants cannot be neglected, even if the final 3D-1D abun- 
dance correction, A3D-1D, is seemingly very small. 

The mean 3D LTE barium-to-iron abundance ratio obtained 
for the eight red giants is ([Ba/Fe])3D lte = 0.28 + 0.07 + 
0.10. The 3D LTE barium abundance measurements made for 
a given star from the three barium lines always agree to within 
* 0.03 dex. In the case of all twenty giants studied here, the 
mean ID NLTE barium-to-iron ratio corrected for the 3D-rela ted 
effects is ([Ba/Fe])3D+NLTE = 0.10 + 0.08 + 0.10 and there- 
fore is only slightly different from the ID NLTE value obtained 
in Sect. 13.31 However, the positive sign of the 3D-1D abun- 
dance differences indicates that in the spectra of red giants in 
NGC 6752 the three studied Ba II lines will be weaker in 3D 
than in ID, in cont rast to wh at is generall y seen in re d giants at 
this metallicity (cf. lCoUet e? al. 2007; Dobrovolskas et al. 2010). 

For the Ba II lines, the 3D-1D abundance corrections are 
sensitive to the choice of microturbulence velocities in the 
ID models: an increase in the microturbulence velocity by 
0.10km s ' leads to an increase of 0.07 dex in the 3D-1D abun- 
dance correction. At the same time, the ID abundance itself de- 
creases by roughly the same amount. The result is that although 
the 3D correction is sensitive to microturbulence, the 3D cor- 
rected abundance is not. 

5. Conclusions 

We have derived the ID LTE and ID NLTE abundances of bar- 
ium for 20 red giant stars in the globular cluster NGC 6752. The 
mean barium-to-iron abundance ratios are ([Ba/Fe])iD lte = 
0.24 + 0.05 + 0.08 and <[Ba/Fe])iD nlte = 0.05 + 0.06 + 0.08 
(the first error measures the star-to-star variation in the abun- 
dance ratio and the second is the systematic uncertainty in the 
atmospheric parameter determination, see Sect. 13.1) . Individual 
barium-to-iron abundance ratios show little star-to-star variation, 
which leads us to conclude that there is no intrinsic barium abun- 
dance spread in the RGB stars at or slightly below the RGB 
bump in NGC 6752. This conclusion is in line with the results 
obtained in other studies, fo r stars in both this and othe r GGCs 
dNorris & Da Cost a 1995; J ames et al.ll2004llYong et ag r2005). 

The derived ID NLTE barium-to-iron abundance ratio is 
comparable t o the one observed in Galactic halo stars of the same 
metallicity ( Sneden et al.l 12008"). It is also similar to the mean 
barium -to-iron abundance ratio obtained by Mish enina et al.l 
( I2OO9I) for 2 red giants in the Galactic globular cluster MIO. We 
therefore conclude that the barium-to-iron abundance ratios ob- 
tained here generally agree with those seen in the oldest Galactic 
populations and are not very different from those observed in 
halo stars. 

We have also obtained 3D LTE barium abundances for 8 red 
giants on the lower RGB in NGC 6752. The mean 3D LTE bar- 
ium abundance, ([Ba/Fe])3D lte - 0.28 + 0.07 + 0.10, is only 
0.05 dex higher than that obtained for these stars in ID LTE. 
This small 3D-1D correction leads to very minor adjustment of 
the mean ID NLTE barium-to-iron ratio for the 20 investigated 
giants, ([Ba/Fe])3D+NLTE = 0.10 + 0.08 + 0.10. 

It would be misleading, however, to conclude that the role of 
the 3D effects in the formation of the barium lines in the atmo- 
spheres of red giants in NGC 6752 is minor As a matter of fact, 
we have found that the 3D-1D abundance corrections owing to 
horizontal temperature inhomogeneities in the 3D model (i.e., 
A3D-(3D> correction) and differences in the temperature profiles 
between the average (3D) and ID models (A(3d>-id correction) 
are substantial and may reach ~ +0.1 dex (Table |5]l. However, 



their sign depends on the line strength, and owing to this sub- 
tle fine-tuning their sum is significantly smaller, from -0.03 to 
0.02 dex, which for this given set of atmospheric and atomic line 
parameters maintains the size of the 3D-1D abundance correc- 
tions at the level of the errors in the abundance determination. 
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Appendix A: Spectral lines Table A.l. List of iron lines used for atmospheric parameter de- 

termination of the target stars. 



Species 


A, nm 


;r,ev 


log«/" 


Fel 


585.22187 


4.548 


-1.330 


Fel 


586.23530 


4.549 


-0.058 


Fel 


590.56720 


4.652 


-0.730 


Fel 


591.62474 


2.453 


-2.994 


Fel 


592.77891 


4.652 


-1.090 


Fel 


593.01799 


4.652 


-0.230 


Fel 


593.46549 


3.928 


-1.170 


Fel 


595.27184 


3.984 


-1.440 


Fel 


597.67750 


3.943 


-1.310 


Fel 


602.40580 


4.548 


-0.120 


Fel 


602.70509 


4.076 


-1.089 


Fel 


605.60047 


4.733 


-0.460 


Fel 


606.54822 


2.608 


-1.530 


Fel 


607.84910 


4.795 


-0.424 


Fel 


608.27106 


2.223 


-3.573 


Fel 


609.66653 


3.984 


-1.930 


Fel 


612.79066 


4.143 


-1.399 


Fel 


613.66153 


2.453 


-1.400 


Fel 


613.69947 


2.198 


-2.950 


Fel 


613.76917 


2.588 


-1.403 


Fel 


615.16181 


2.176 


-3.299 


Fel 


615.77284 


4.076 


-1.260 


Fel 


617.33356 


2.223 


-2.880 


Fel 


618.02042 


2.727 


-2.586 


Fel 


618.79904 


3.943 


-1.720 


Fel 


619.15584 


2.433 


-1.417 


Fel 


620.03129 


2.608 


-2.437 


Fel 


621.34303 


2.223 


-2.482 


Fel 


621.92810 


2.198 


-2.433 


Fel 


623.07230 


2.559 


-1.281 


Fel 


623.26412 


3.654 


-1.223 


Fel 


624.06462 


2.223 


-3.233 


Fel 


624.63188 


3.602 


-0.733 


Fel 


625.25554 


2.404 


-1.687 


Fel 


626.51340 


2.176 


-2.550 


Fel 


627.02250 


2.858 


-2.464 


Fel 


630.15012 


3.654 


-0.718 


Fel 


632.26855 


2.588 


-2.426 


Fel 


633.53308 


2.198 


-2.177 


Fel 


633.68243 


3.686 


-0.856 


Fel 


634.41491 


2.433 


-2.923 


Fel 


635.50290 


2.845 


-2.350 


Fel 


638.07433 


4.186 


-1.376 


Fel 


639.36013 


2.433 


-1.432 


Fel 


640.00012 


3.602 


-0.290 


Fel 


641.16493 


3.654 


-0.595 


Fel 


641.99496 


4.733 


-0.240 


Fel 


642.13508 


2.279 


-2.027 


Fel 


643.08464 


2.176 


-2.006 


Fel 


647.56244 


2.559 


-2.942 


Fel 


648.18703 


2.279 


-2.984 


Fel 


649.49805 


2.404 


-1.273 


Fel 


649.64666 


4.795 


-0.570 


Fel 


651.83671 


2.831 


-2.460 


Fel 


659.38705 


2.433 


-2.422 


Fel 


660.91103 


2.559 


-2.692 


Fel 


663.37497 


4.558 


-0.799 


Fel 


667.79870 


2.692 


-1.418 


Fel 


670.35674 


2.758 


-3.160 


Fel 


675.01525 


2.424 


-2.621 


Fel 


680.68449 


2.727 


-3.210 


Fell 


599.13760 


3.153 


-3.540 


Fell 


608.41110 


3.199 


-3.780 


Fell 


614.92580 


3.889 


-2.720 


Fell 


624.75570 


3.892 


-2.310 


Fell 


636.94620 


2.891 


-4.160 


Fell 


641.69190 


3.892 


-2.650 


Fell 


643.26800 


2.891 


-3.520 


Fell 


645.63830 


3.903 


-2.100 


Fell 


651.60800 


2.891 


-3.320 



